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Abstract 

Conventional weak-coupling perturbation theory suffers from problems that 
arise from resonant coupling of successive orders in the perturbation series. 
Multiple-scale perturbation theory avoids such problems by implicitly per- 
forming an infinite reordering and resummation of the conventional perturba- 
tion series. Multiple-scale analysis provides a good description of the classical 
anharmonic oscillator. Here, it is extended to study the Heisenberg opera- 
tor equations of motion for the quantum anharmonic oscillator. The analysis 
yields a system of nonlinear operator differential equations, which is solved 

exactly. The solution provides an operator mass renormalization of the theory. 
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Multiple-scale perturbation theory (MSPT) is a powerful and sophisticated perturbative 
method for solving physical problems having a small parameter e ]1]J2|]. MSPT is applicable 
to both linear and nonlinear problems. Indeed, it is so general that other perturbative 
methods such as WKB theory and boundary-layer theory, which are useful in more limited 
contexts, may be viewed as special cases of MSPT JT| . 

MSPT recognizes that dynamical systems exhibit characteristic physical behaviors at 
various length or time scales. The problem with conventional perturbation theory is that 
there is often a resonant coupling between successive orders. This coupling gives rise to 
secular terms (terms that grow rapidly with the length or time variable) in the perturbation 
series. Secular terms conflict with physical requirements that the solution be finite. MSPT 
reorganizes the conventional perturbation series to eliminate secular terms, and in doing so 
it describes quantitatively the behaviors that occur at many scales. 

Ordinarily, MSPT is applied to classical differential equations such as Duffing's equation 
(the nonlinear equation of motion for the classical anharmonic oscillator) : 

y" + y + 4ey 3 = (e>0). (1) 

The positivity of e ensures that y(t) is bounded JTJ. The classical harmonic oscillator (e = 0) 
has only one time scale, the period of oscillation. However, when e 7^ 0, the nonlinear term 
in Eq. ([!]) introduces many time scales. Using MSPT one can show that on a long-time scale 
[t = Oie- 1 )] there is a frequency shift of order e. 

In this letter we generalize MSPT and apply it to the Heisenberg operator equations of 
motion of the quantum anharmonic oscillator (the quantum version of Duffing's equation). 
This generalization of MSPT techniques is nontrivial because it gives rise to a nonlinear 
system of operator differential equations || . We find the exact closed-form solution to this 
system and thereby obtain the quantum operator analog of the classical frequency shift - 
an operator mass renormalization that expresses the first-order shift of all energy levels. 

We begin our presentation by reviewing the difficulties one encounters when one tries to 
solve Duffing's equation using a conventional perturbation series, y(t) = J2™=o e "?/n(i)- We 
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choose as initial conditions 



y(0) = l and y'(0) = 0, (2) 

which translate into y n (0) = ^n,o an d 2/n(0) = 0' an d substitute into Eq. ([!]). To zeroth- 
and first-order in powers of e, we have 

2/o+2/o = 0, (3) 
y'; + yi = -Ay 3 . (4) 

The solution to Eq. (§) satisfying the initial conditions is yo(t) = cost. Introducing this 
solution into Eq. (f|), we obtain y" + yi — — cos(3t) — 3 cost, which represents a forced 
harmonic oscillator whose driving term has frequencies 3 and 1. A harmonic oscillator when 
driven at its natural frequency, which in this case is 1, exhibits resonance. As a result, 
the solution yi(t) = | cos(3t) — |cost — |tsint contains a secular term that grows linearly 
with time t. The function yi(t) cannot be valid for long times because the exact solution 
to Dufnng's equation remains bounded for all t jl]. Hence, the conventional perturbation 
expansion is sensible only for short times t << e -1 . 

If one is clever, one can use the conventional perturbation series to determine y(t) for 
long times, say of order e~ 1 . To do so we note |l| that the structure of the most secular 
(highest power in t) term in y n {t) has the form ^(3it/2) n e tt /n\ + c.c; we then approximate 
y(t) by summing the most secular term in every order and the result is a cosine function 
that remains bounded for all times t: 
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Hence, on the long-time scale r = et we see a frequency shift in the oscillator of order |e. 
This result is not exact because there are less secular terms to all orders in the perturbation 
expansion; such terms give rise to frequency shifts of order e 2 , e 3 , ... . 

The advantage of MSPT is that it reproduces Eq. (|5]) directly and bypasses the elaborate 
procedure of summing the conventional perturbation series to all orders by excluding ab initio 



secular terms from the perturbation expansion. MSPT assumes the existence of many time 
scales (t, t = et, a = e 2 t, ...), which can be temporarily treated as independent variables. 
Here, we illustrate by performing just a first-order calculation. We use only the two variables 
t and t = et and seek a perturbative solution to Eq. (|1|) of the form 

y(t)=Y (t,T)+eY 1 (t,r) + O(e 2 ). (6) 

The chain rule and the identity % = e convert Eq. (JED to a sequence of partial differential 
equations for the dependent variables Y , Y\, ... . The first two are 

l^o + Y« = 0. (7) 

^y 1 + n = -4y3-2^y„. (8) 

The general solution to Eq. (0) is Y (t, r) = A(t) cost + B(t) sint. We substitute Y (t, r) 
into the right side of Eq. (j§) and use triple-angle formulas such as cos 3 1 — \ cos(3t) + | cost 
to simplify the result. To determine the functions A(j) and B(r) we demand that there 
be no resonant coupling between zeroth and first order in perturbation theory so that no 
secular terms appear in Y\. That is, we require that the coefficients of sint and cost vanish: 

2^ = -3A 3 - 3AB 2 and 2^ = 3B 3 + 3A 2 B. (9) 
dr dr 

To solve this system we multiply the first equation by B(t), the second by A(t), and 
add the resulting equations. Letting C(r) = ~ [^4(t)] 2 + \ [B(t)} 2 , we obtain 

^ W = o. (10) 

Thus, C(t) is the constant C(0) and the differential equation system @ becomes linear: 

^-B = -3C(0)A and ^—A = 3C(0)B. (11) 
dr dr 

The solution to this system that satisfies the initial conditions is C(0) = \ and Yo(t,T) = 



cos 



[1 + *e)t 



which reproduces the approximate solution in Eq.(^). While conventional 



perturbation theory is valid for t « e 1 , the MSPT approximation is valid for t « e 2 . 



We now consider the quantum anharmonic oscillator, whose Hamiltonian is H(p, q) = 
\v 2 + \<f + e <? 4 - Here, e > so that the spectrum of H (p, q) is bounded below and p and q 
are operators satisfying the canonical equal-time commutation relation [q(t),p(t)] = ih. 

The Heisenberg operator equations of motion, ^q = ■ifg, H(p,q)] = p and -^p = 
■4[p, H(p, q)] = — q — 4eq 3 , combine to give the quantum Duffing's equation (|T|): 

d 2 

—q + q + 4eq 3 = 0. (12) 

Since p(t) and q(t) are operators, we cannot impose numerical initial conditions like those 
in Eq. (0); rather, we enforce a general operator initial condition at t = 0: 

q(0)=q and p(0) = p Q , (13) 

where p and q are time-independent operators obeying the Heisenberg algebra [(ftbPo] — ih- 
We now apply MSPT to Eq. ([T2]). Assuming that q(t) exhibits characteristic behavior 
on the short-time scale t and on the long-time scale r = et, we write 

q(t) = Q(t, r) = Q (t, r) + eQ^t, r) + 0(e 2 ). (14) 

This equation is analogous to Eq. (^) but here Qo and Q\ are operator-valued functions. 

We substitute q{t) in Eq. ( ITj ) into Eq. flT2"D, collect the coefficients of e° and e 1 , and 
obtain operator differential equations analogous to Eqs. (0) and (H): 

^Qo + Qo = 0, (15) 
^Qi + gi = -4Q^-2^-g . (16) 

Because Eq. (|l|) is linear, it is easy to find its general solution, 

Qo(t,T) =A(t) cost + B(t) sin t, (17) 

and from p = % we obtain the momentum operator p(t) = B(r) cost — A(t) sint + 0(e). 

It is now necessary to find the coefficient functions A(t) and B(r), which are operators. 
The canonical commutation relation in [q(t),p(t)] = ih implies that these operators satisfy 
[A(t), B(t)} = ih. Also, the initial conditions in Eq. (|13|) give 



A(0) = go and 5(0) = p - (18) 

We determine A(t) and B{t) by evaluating the right side of Eq. (p!|) and expanding the 
cubic term, taking care to preserve the order of operator multiplication. As in the classical 
case, we simplify the result using triple-angle formulas. To prevent secularity in Qi(t, r) we 
set the coefficients of cos t and sin t to zero and obtain 

2— = -3„4 3 - BAB - BBA - ABB, 

dr 

d A 

2-— = 3£> 3 + ABA + AAB + BAA. (19) 
dr 

This system of operator- valued differential equations is the quantum analog of Eq. (|9|). 

To solve the system fllTf ) we pre- and post-multiply the first equation by B{j) and the 
second equation by A{j). Adding the resulting four equations, we get 

= < 20) 

where Ti = \A 2 + \B 2 . Equation (|2"0|), the quantum analog of Eq. (|T0"D, shows that 7i is 
independent of the long-time variable r. Thus, Eq. ([T^) allows us to express TL in terms 
of the fundamental operators po and qo: H = \p\ + \q\. We then use the commutator 
\A(t\B(t)\ = ih to rewrite Eq. (|19|) in manifestly Hermitian form: 

JL B = -1(HA + AH), ^A=l(HB + BH). (21) 
dr l dr A 

Suppose for a moment that we could replace the operator 7i by the numerical constant 

C(0) in Eq. (0). Then we would obtain the c-number coupled differential equations in 

Eq. (|TTD. That system is linear so we could ignore operator ordering and easily find the 

solution satisfying the initial conditions ([18]): A{r) = q cos[3C(0)r] + p sin[3C(0)r] and 

B{t) = po cos[3C(0)t] — g sin[3C(0)r]. This solution suggests the structure of the exact 

solution to the operator differential equation system (pl|). The formal solution is a natural 

generalization using Weyl-ordered products of operators: 

A(t) = W[g cos(3ftr) + p sin(3Wr)], 

B(t) = W[p cos(3Ht) - g sin(3Wr)]. (22) 



The operator ordering W[<7q/(Wt)] is defined as follows: (1) Expand f(Hr) as a Tay- 
lor series in powers of the operator 7ir; (2) Weyl-order the Taylor series term-by-term: 
W{q H n ) = ^E?=o(i)^ i <?oH n -' ? '. Using this definition it is easy to verify that Eq. @ 
is indeed the exact operator solution to Eq. ( BTj ) satisfying the initial conditions (0). We 
simplify the formal solution in Eq. fl22| ) and re-express it in closed form by observing that if 
we reorder W(qoli. n ) by commuting q symmetrically to the left and to the right to maintain 
the Hermitian form we generate a set of polynomials @] of degree n: 

h n 



W(q H n ) o 



We identify E n as the nth Euler polynomial || in which the argument is shifted by \: 



H H 2 _ 1 7Y 3 3 H 3H 2 5 



and so on. The generating function for these nonorthogonal polynomials is 

^Ti-^S^U + 2-) (W<')- 

This generating function allows us to express the following Weyl-ordered product compactly: 

Using Eq. (^3j) we rewrite compactly the cosines and sines in Eq. (0), substitute this 
result into Eq. fllTD , replace r by et, and obtain 

q cos(t + 3Het) + cos(t + 3Het)q Q 



Qo(t,r) 



2cos(3etft/2) 
p sin(£ + 3TCet) + sin(t + 3Het)p 



(24) 



2cos(3e^/2) 

Equation (j24"D is the quantum operator analog of Eq. (H) and is the objective of our multiscale 
analysis. [We can recover the classical MSPT approximation in Eq. fl5|) by taking the limit 
h —>■ and imposing the classical initial conditions po = and go — lj which give 7i — |.] 
Recall that in the classical case we identify the coefficient of the time t as a first-order 
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approximation to the frequency shift. Since the coefficient of t in Eq. fl24|) is an operator, 
we have derived an operator form of mass renormalization. 

To elucidate this operator mass renormalization we study matrix elements of Eq. (pl|). 
The time dependence of a matrix element reveals the energy-level differences of the quantum 
system. It is easy to construct a set of states because the operators q and p satisfy the 
Heisenberg algebra [go,.Po] — ih- Hence, appropriate linear combinations of go an d Po may 
be used as raising and lowering operators to generate a Fock space consisting of the states 
|n). By construction, these states are eigenstates of the operator H: TC\n) = (n + |) h\n). 
Evaluating Eq. (|2~4D between the states (n — 1| and \n) and allowing the operator 7i to act 
to the left and the right, we obtain 

(n — 1|Q |^) — (n — l\qo\n) cos[t(l + 3nhe)} 

+(n-l\po\n)sm[t(l + 3nhe)], (25) 

which predicts that the energy-level differences of the quantum oscillator are 1 + 3nhe. We 
may verify this result by recalling that the first-order correction to the energy eigenvalues 
@ is E n = n + | + \eh(2n 2 + 2n + 1) + 0(e 2 ). Thus, E n - £ n -i = 1 + Me + 0(e 2 ). 

We conclude by noting that in addition to the operators q(t) and pit), the wave function 
i/j(x) for the quantum anharmonic oscillator also exhibits multiscale behavior. Specifically, 
to all orders in conventional weak-coupling Rayleigh-Schrodinger perturbation theory, ijj(x) 
behaves like the Gaussian e _x ' 2 / 4 for large x; however, a geometrical-optics approximation 
to ip(x) from WKB theory predicts that for large x, ip(x) decays like the exponential of a 
cubic e-^N'/e. We 

can resolve this discrepancy by using a MSPT approach; we reorder the 
perturbation series by resumming secular terms 0. 

The wave function ip(x) obeys the Schrodinger equation 

(-^ + \x 2 + \^ - E{e)\ j iP(x) = (26) 

and satisfies ^(±oo) = 0. The conventional perturbative approach to Eq. (^) [^] represents 
both the eigenfunction and eigenvalue as asymptotic series in e: ip(x) ~ S^Lo e ""yn(x) and 
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E(e) ~ Y^=o£ n E n - In Ref. it is shown that for the ground state, y n (x) is a Gaussian 
multiplied by a polynomial of degree 2n in the variable x 2 , y n = e~ x2 ^ A P n (x), where 

2n i 

P {x) = l and P n ( x ) = Y J C n , k (--x 2 ) k (n > 0). 

i=i ^ 

The recursion relation for the polynomials P n (x) is typical of all perturbative calculations; 
the homogeneous part of this recursion relation is independent of n while the inhomogeneous 
part contains all previous polynomials. This recursive structure is responsible for succes- 
sive orders of perturbation theory being resonantly coupled and causes the degree of the 
polynomials to grow with n. 

To study the behavior of the wave function ip(x) for large x we approximate i/j(x) by 
resumming the perturbation series and keeping just the highest power of x in every order. 
[This is an exact analog of finding the coefficient of the highest power of t (most secular 
term) in nth order in perturbation theory for the classical anharmonic oscillator.] This 
resummation gives a new representation of ip(x): e - x2 /^ e - ex4 / 16 multiplied by a new set of 
polynomials. For the classical anharmonic oscillator summing leading secular terms also gives 
an exponential approximation [see Eq. (|5|)]. However, the classical and quantum anharmonic 
oscillators are quite different; although we have summed the most secular terms to all orders 
in perturbation theory, the result is not the actual behavior of the wave function ip(x) for 
large x\ the correct behavior is an exponential of a cubic and not a quartic! 

If we iterate this resummation process we find that each reorganization of the perturba- 
tion series gives an additional term in the exponential 

exp _ + _L f v - J^v + ...) (27, 

and a new set of polynomials. It seems impossible for this approach to give cubic exponential 
behavior because at each stage in the reorganization of the perturbation series the variable 
x appears only in even powers. However, we recognize that the exponent in Eq. ( p?|) is the 
beginning of a binomial series whose sum is ^ 1 — (1 + ex 2 ) 3 / 2 . If we now let x be large 
(ex 2 » 1), we recover the correct asymptotic behavior of ip(x) ||. 

9 



The approach used above for the anharmonic oscillator wave function has been used 
in perturbative quantum field theory to sum leading-logarithm divergences || and lead- 



ing infrared divergences [|Kj. It is our hope that in the future the direct nonperturbative 
multivariate approach of MSPT will provide a framework to simplify such schemes. 
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